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Summary. The basis for a hydrodynamic description of granular gases is discussed 
for a low density gas of smooth, inelastic hard spheres. The more fundamental meso- 
scopic description is taken to be the nonlinear Boltzmann kinetic equation. Two ar- 
guments are presented in favor of a hydrodynamics for granular gases. The first one is 
the concept of a "normal" solution and its explicit approximate construction via the 
Chapman-Enskog method. The second is the demonstration of hydrodynamic modes 
in the spectrum of the generator for the dynamics of small spatial perturbations of 
the homogeneous reference state. In the first derivation of the nonlinear hy- 

drodynamic equations is given to Navier- Stokes order, with explicit expressions for 
the transport coefficients. The approach is formal and the context of the derivation 
is left unspecified, although internal mathematically consistency is demonstrated. 
The second method is more restricted, leading only to linearized hydrodynamics, 
but with the potential to define more sharply the context of hydrodynamics. 



1 Introduction 

The form and validity of a hydrodynamic description for granular gases re- 
mains a controversial subject [1]. The objective here is to formulate the con- 
ceptual and mathematical questions in terms of the underlying Boltzmann 
kinetic equation, and to describe some recent results providing support for 
hydrodynamics [2, 3, 4]. Two distinct approaches will be discussed. The first 
is a consideration of "constitutive equations" for a closure of the exact macro- 
scopic balance equations for the average mass, energy (or temperature), and 
momentum (or flow velocity). These equations follow directly from the exis- 
tence of a "normal" solution to the Boltzmann equation. The expected con- 
ditions for a normal solution, e.g., away from spatial or temporal boundary 
layers, characterize the context in which a hydrodynamic description can be 
expected. The explicit construction of an approximate normal solution under 
conditions of small spatial variation of the hydrodynamic fields is given by 
an extension of the Chapman-Enskog method [5]. It is emphasized that no 



2 



James W. Dufty and J. Javier Brey 



additional assumptions beyond the normal form and small gradients are re- 
quired by the inelasticity of collisions. Differences from normal gases (e.g., the 
reference state, linear stability, form of the constitutive equations) are a con- 
sequence of the method, not a limitation on it. The resulting hydrodynamics 
comprises the Navier- Stokes approximation for a granular gas. Finally, it is 
emphasized that a breakdown of the latter does not imply the absence of a hy- 
drodynamic description. Indeed, it is an interesting new challenge of granular 
gases to obtain constitutive equations for more general conditions [1, 6]. 

The second means to explore the existence and context of hydrodynamics 
is to consider small spatial perturbations of a strictly homogeneous state. This 
problem is somewhat simpler as it is governed by the linearized Boltzmann 
equation, and it is then sufficient to consider a single Fourier component. The 
spectrum of the generator for the dynamics determines all possible excitations, 
and it can be asked if the hydrodynamic excitations are among these. For 
consistency, the hydrodynamic excitations found in this way must be the 
same as those from the Navier-Stokes equations under the same conditions. 

Recently, the eigenvalue problem for this generator has been considered 
[2, 3, 4]. There is a technical problem associated with the "cooling" of the 
reference homogeneous state, which is controlled by the introduction of suit- 
able dimcnsionlcss variables. This entails a corresponding change in the form 
of the generator for the dynamics. The existence of hydrodynamic modes has 
been proved by the exact construction of the hydrodynamic eigenvalues and 
eigenfunctions in the long wavelength limit. Their form to Navier-Stokes or- 
der is then obtained by perturbation theory, and complete agreement with the 
results from the Chapman-Enskog method is found. The dominance of these 
hydrodynamic excitations at long times requires that their eigenvalues be the 
smallest ones and be isolated from the rest of the spectrum. Some comments 
are offered on why this is to be expected. Elsewhere in this volume [7] these 
expectations are confirmed for an exactly solvable kinetic model. 

The objective of this brief overview is to clarify the concept of a hydro- 
dynamic description in its most general sense, and to sharpen the sufficient 
conditions for its validity. 

2 Boltzmann Equation and Balance Equations 

Consider a system of N smooth hard spheres of mass m and diameter a in a 
large volume V. If the density n = N/V is sufficiently small, na 3 << 1, the 
phase space density, /(r, v, t), for the number of particles with position r and 
velocity v at time t is determined from the Boltzmann equation [8, 9, 10], 



where J is the inelastic Boltzmann collision operator. The most important 
properties following from this collision operator for the following discussion 
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are: 

1 \ / 



/ 



dv| v J[r,v|/(t)]= |. (2) 

i m (v-u) 2 / \-%n(r,t)T( T ,t)C(T,t) , 

The two zeros on the right side of (2) follow from conservation of mass and 
momentum. The last term results from non-conservation of energy, where £ 
is the fractional rate of decrease in the energy due to the inelasticity, 

«'■ " = 2MM^) / * / )«'• • «>■ < 3 » 

Here g = v — vi and a is the restitution coefficient (a = 1 is the elastic limit). 
The "hydrodynamic fields" appearing in (2) and (3) are the density n(r,t), 
the temperature T(r,t), and the macroscopic flow velocity u(r, t). These are 
defined in terms of velocity moments of /(r, v, t) by 

n(r,t) \ . ( 1 \ 
n(r,t)u(v,t) = / dv v /(r,v,t). (4) 

§n(r,t)T(r,t)/ 7 ^ I™ (v - u) 2 / 

The macroscopic balance equations for these fields follow directly from the 
Boltzmann equation using the properties (2), 

d t n + V • (mi) = 0, (5) 

(& + u-V)ui + (mn) -1 Vji^=0, (6) 

(9 t + u • V + C) T + J- (Py Vj-Ui + V • q) = 0. (7) 

The pressure tensor and heat flux q are defined by 

Pij(r,t) = J dvm(vi - Ui)(vj - Uj)f(r,v,t), (8) 

dv-(v-u) 2 (v-u)/(r,v,i). (9) 

In the above equations and in the following, summation over repeated indexes 
is implicit. For a normal gas with elastic collisions, Eqs. (5)-(9) are the local 
macroscopic conservation laws for mass, momentum, and energy. The only for- 
mal difference here is the appearance of the "cooling rate", £, in the equation 
for the temperature. Moreover, all the explicit dependence on the restitution 
coefficient a occurs in this cooling rate. 

The macroscopic balance equations are the starting point for developing 
a hydrodynamic description of the gas. They are an exact consequence of the 
Boltzmann equation, but are of limited utility without further specification 
of the three unknown functions P^ , q, and £. Since these are specific func- 
tionals of /(r, v,t), the above requires a solution to the Boltzmann equation. 
The special class of solutions leading to hydrodynamics, normal solutions, arc 
described in the next section. 
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3 Normal Solution, Constitutive Equations, and 
Hydrodynamics 

A sufficient condition for a hydrodynamic description is the existence of a 
normal solution to the Boltzmann equation. To describe it, denote the hy- 
drodynamic fields collectively by the set {y a (r, t)}. Then, a normal solution 
is defined as one for which all space and time dependence of the distribution 
function is determined by the hydrodynamic fields, i.e., 

/(r,v,t) =f[r,v,t\{y a }], (10) 

such that 

If such a solution can be found, its use in the functional (3, (8), and (9), 
converts them to normal forms as well, 

C = C[r,t I {»«}], Pij =P ij [r,t\{y a }}, q = q [r, t | {y a }} . (12) 

Such normal forms are called "constitutive" equations. Finally, use of these 
normal forms in the macroscopic balance equations gives a closed set of equa- 
tions for the hydrodynamic fields: 

dty a (T,t)=N a [r,t\{yp}], (13) 

where N a [r, t | {yp}] is, in general, a non-local, non-linear functional of the 
set of fields {yp}. Equations (13) are the hydrodynamic equations in their 
most general form, encompassing both rheological and viscoelastic nonlin- 
ear transport. Such a general context is typically reserved for complex fluids 
(e.g., colloids, large molecular weight fluids) and is of little interest for normal 
gases. However, structurally simple granular gases can exhibit many of the 
non-Newtonian properties of complex fluids. Therefore, in discussing the ap- 
plicability of a hydrodynamic description for granular gases, it is important to 
admit the possibility of a closed set of equations more complex than the local 
partial differential equations of the Navier-Stokes approximation (see below). 

The construction of a hydrodynamic description is seen to be intimately 
connected to finding a normal solution to the Boltzmann equation. In fact the 
two problems comprise a single self-consistent problem. For a normal solution, 
the Boltzmann equation (1) becomes 

= J[r,v|/], (14) 

For specified fields, this is an equation for the velocity dependence of the nor- 
mal distribution as a functional of the fields. This dependence then allows 
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determination of the normal forms in (12), and hence N a [r,t \ {yp}]- Finally, 
solution of Eq. (13), with suitable initial and boundary conditions, provides 
the explicit forms for the fields, and completes the normal solution. The exis- 
tence and determination of this solution is the central problem for establishing 
a hydrodynamic description for both normal and granular gases. The concept 
of a normal solution and its use in the macroscopic balance equations makes 
no special reference to the possible inelasticity of collisions. 

To clarify the conditions under which a normal solution could be expected, 
consider first a normal gas with elastic collisions in an initial non-equilibrium 
state with specified hydrodynamic fields {y a (r, t = 0)}, whose values vary 
smoothly across the system. In each small region (e.g., cells of 20-30 parti- 
cles), the velocity distribution approaches a Maxwellian characterized by the 
hydrodynamic fields at its central point r. Subsequently, exchange of mass, en- 
ergy, and momentum tends to equilibrate these fields to uniform values (or to 
steady values if the system is driven). The first stage, approach to a universal 
form for the velocity distribution, occurs after a few collisions. This estab- 
lishes the normal form of the solution, where the hydrodynamic fields and 
their gradients characterize the state. Deviations from the Maxwellian are 
due to fluxes of mass, momentum, and energy across the cells. These fluxes 
are proportional to the differences in values of the fields (i.e., to their spatial 
gradients). The second stage is the evolution of the distribution through the 
changing values of the fields, according to the hydrodynamic equations. 

This two-stage evolution can be expected for granular gases as well. The 
initial velocity relaxation will not approach the local Maxwellian, but some 
other corresponding normal distribution determined from the inelastic Boltz- 
mann equation (see below) . Subsequently, the deviations from this normal dis- 
tribution characterizing spatial inhomogeneities, will again be via the macro- 
scopic balance equations for the granular gas. 

3.1 Constructive Example: Navier-Stokes Approximation 

It is easily verified that the normal solution cannot be simply a local function 
of the fields alone, i.e. / ^ /[v | {y a }]- However, the nonlocal functional 
dependence can be converted to a local function of the fields and all their 
gradients, 

/ = / (v, {y a (r, *)} , {Vy a (r, *)} , { VVy Q (r, t)} , . . .) . (15) 

The time derivatives have not been included, since they can be expressed for- 
mally in terms of the gradients through the hydrodynamic equations (13), with 
the corresponding local form N a ({y f3 (r, t)} , {Vyp (r, t)} , { Wyp (r, t)} ,...). 
Such a representation is useful when some higher order gradients can be ne- 
glected. For example, if the initial and boundary conditions suggest nearly 
linear hydrodynamic profiles (such as for shear flow between parallel moving 
plates), the above can be estimated as / ~ / (v, {y a (r, t)} , {Vy a (r, t)}). This 
approximation could still include large gradients and nonlinear transport. 
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A simpler case is that for which all the spatial gradients are small. More 
precisely, consider states for which the fractional changes in the hydrodynamic 
fields are small over a mean free £, i.e., £|Vy Q (r, t) jy a (r, t)\ <~ e << 1 and 
£ | V n+1 y a (r, t) /V n y a (r, t) | ~ e << 1. The parameter e and is called the uni- 
formity parameter. The normal solution can now be calculated perturbatively 
as an expansion in this small parameter. To first order in e, it must be linear 
in the hydrodynamic gradients 

/(r, v, t) ~ /(°) (v, {y a (r, t)} ) (l + eX^ [v, {y, (r, i)}] • (r, t) ) . (16) 

Before discussing the specific forms for /(°) and X Q , the normal forms for 
£(r, t, ), Pij(r,t,), and q(r, i, ) can be determined directly by use of (16) in 
the functional (3), (8), and (9), and invoking the fluid symmetry, 

C=(l-a 2 )n(r,t)a 2 (^My /2 C, (17) 
q = -k [T (r, t)] VT (r, t) — fi[T (r, t) , n (r, t)] Vn (r, t) , (18) 



P ij =p(v,t) 6ij-ri[T(T,t)] 



2 

-%V-u(r, t) 



Vjtij (r,t) + V»tij (r,t) 



(19) 



Here p = nT is the pressure, £ is a dimensionless constant determined from 
/(°), and the transport coefficients k, fi, and r\ are linear functionals of the X Q . 
Their exact expressions in terms of Z^ ' and the X a are given elsewhere [5, 11] 
and will not be repeated here. In the expression for the cooling rate (17), only 
the leading contribution (zeroth order in the gradients) is given. In principle, 
consistency would require to consider up to second order in the gradients, since 
this is the order of the terms associated with the heat flux and the pressure 
tensor. Nevertheless, contributions from the cooling rate of first order in the 
gradients of the fields are seen to vanish because of symmetry reasons, and 
those of second order have been shown to be negligible as compared with the 
corresponding contributions from Eqs. 18) and (19) [5]. 

The macroscopic balance equations of the last section together with (17)- 
(19) constitute the Navier-Stokes hydrodynamic equations. Their form is the 
same as for a normal fluid, except for the presence of the cooling rate ( and 
the new transport coefficient /i in the equation for the heat flux. Of course, 
the values of the transport coefficients are different, depending on the value 
of the restitution coefficient a. 

The leading term /W in (16) is determined from the Boltzmann equation 
to order e°, 

7 fl/(0) 



s dT 



v|/W . (20) 
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Since /(°) is normal, it follows from dimensional analysis that its dependence 
on temperature can occur only through normalization and a scaling of the 
velocity. Then (20) becomes: 



2 <9v 



(v-u)/W 



J 



r, v 



1/(0) 



(21) 



This equation also arises when considering the Boltzmann equation for a spa- 
tially homogeneous state, the so-called homogeneous cooling state (HCS), as 
discussed in the next section. Here, it is determining the same velocity de- 
pendence, but scaled with respect to the exact hydrodynamic fields for the 
nonuniform, nonstationary state. Thus the solution here is more appropriately 
named the local HCS. It is important to note that this leading order "refer- 
ence" state is determined by the Chapman-Enskog procedure itself, not by an 
arbitrary, independent choice. 

Similarly, the first corrections, determined by X a , are obtained from the 
Boltzmann equation to order e, with the result: 



01n/(°) 



Vy a - X T ■ V (CO 



dy 



7V« (r,t | {yp})+v-Vy a 



(22) 



where Na^ is the term of order e in the expansion of N a , Xt is the coefficient 
of VT in Eq. (16), and L is the linearized Boltzmann collision operator, 



LY(v) 



-/(0)-i (v) / 



dV 



5 J [r,v|/(°)] 



/(°V)r(v'), (23) 



for arbitrary Y(v). By equating coefficients of the hydrodynamic gradients 
in (22), the coefficients having independent scalar, vector, and tensor trans- 
formation properties, a set of linear inhomogeneous integral equations are 
obtained. The Fredholm conditions have been proved for these equations, so 
that existence of solutions is assured. The details are described in [5] and [11]. 
In summary, a normal solution to the Boltzmann equation up through first 
order in e is given by (16). This solution determines the constitutive equations 
in normal form also to this order. These constitutive equations together with 
the macroscopic balance equations, give approximate hydrodynamic equations 
which are the Navier-Stokes equations for a granular gas. The transport co- 
efficients in these equations are also provided by the derivation. 

The main point of the discussion in this section has been to show that 
hydrodynamic equations for granular gases have the same conceptual basis as 
for normal gases: 1) exact macroscopic balance equations for n, T, and u, with 
exact expressions for the unknown fluxes (and the cooling rate) as functionals 
of the distribution function /; 2) a rapid relaxation of / to a normal solution. 
These two ingredients are sufficient for the existence of a hydrodynamic de- 
scription (closed, deterministic equations for the hydrodynamic fields). The 
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construction of the normal solution is a difficult, unsolved problem in gen- 
eral. However, for states with small spatial gradients, the Chapman-Enskog 
perturbation expansion leads to an explicit approximate normal solution. No 
explicit reference to inelasticity occurs at this level. The detailed computation 
of the normal solution and characterization of the space and time scales for 
its validity does, of course, entail differences from that for normal gases. But 
these are technical issues that should be separated from the conceptual basis. 

Property 2) above is the most uncertain component of the discussion. If a 
normal solution exists, what is its context? For molecular gases, the assump- 
tion is that it applies only outside certain "boundary layers": a few mean 
free paths away from spatial boundaries, a few mean free times later than 
any imposed initial conditions, and away from "discontinuities" (e.g., shock 
fronts). Similar restrictions are expected to apply for granular gases, and it 
may be that these are stronger constraints in that case. In addition, there 
is a new time scale for granular gases set by the cooling rate. The hydro- 
dynamic equations described here include the dynamics of cooling, so there 
are new questions regarding the separation of hydrodynamic from all other 
time scales. Monte Carlo and molecular dynamics simulations suggest that the 
separation remains valid even at strong dissipation, but it remains to demon- 
strate this directly from the Boltzmann equation. The next section provides 
a more restricted context in which this question can be posed clearly. 



4 Linearized Boltzmann Equation and Hydrodynamic 
Modes 

In this section, the existence of hydrodynamics for granular gases is considered 
from a different perspective, independent of the concept of a normal solution 
to the Boltzmann equation. The idea is to study the initial value problem 
for small perturbations of the homogeneous state of an isolated system. In 
this case, the Boltzmann equation can be linearized about the homogeneous 
reference state. For appropriate dimensionless variables, the general solution 
is characterized by a linear eigenvalue problem. The eigenvalues determine 
all modes of excitation in the gas due to the small perturbations. If hydro- 
dynamics exists, the hydrodynamic modes must be among these excitations. 
These modes are defined in the long wavelength limit as being the same as 
those obtained from the exact linearized macroscopic equations The reference 
homogeneous state for an isolated granular gas is taken of the form 

fhc (v, t) = nv 3 (t)I* hcs [v/v (t)} , (24) 

with 

Use of this scaling form in the expression for the cooling rate (3), gives again 
the form (17), now as an exact result, 
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at) 



(l-a) 2 ira 2 n (2T{t)\ 1/2 



12 

x 



J dv J (Miv-v.j 3 /;„(»)/;„(»,). (26) 



The macroscopic balance equations for the homogeneous state reduce to 
dtn = dtu = 0, and 

cUnT(i) = -C(*). (27) 

Since the dependence of ((t) on the temperature is known, the time depen- 
dence of the temperature can be calculated. Substitution of (24) into the 
Boltzmann equation gives the equation for fhcs- It is found to be the same 
as that for the local HCS distribution /(°), given by (20) or (21) above. The 
difference here is that the hydrodynamic fields parameterizing this distribu- 
tion are those for the homogeneous state, whereas in the Chapman-Enskog 
context they are the fields for a general inhomogeneous state. As already men- 
tioned, the distribution fhcs is referred to as the HCS solution. Its velocity 
dependence is the same as that of the local HCS solution. 

In the following, it will be assumed that fhcs is known. Consider an initial 
small spatial perturbation of the HCS (other examples might arise from small 
external forces or boundary conditions). Since the perturbation is small ini- 
tially, it is assumed that the deviation of the distribution function from fh cs 
remains small over the times of interest, such that the Boltzmann equation 
can be linearized. Consequently, it is sufficient to consider a single Fourier 
component, 

/ (r, v, t) = fhcs (v, t) [1 + e lkr 5 (k, v, *)] . (28) 

Substitution into the Boltzmann equation and retaining only linear order in 
g gives the linear Boltzmann equation, 



dt+i k-v+L+f- c 1 s (v)^-(vf hcs ) 



ff (k,v,t)=0. (29) 



Here, L is the linear operator defined in (23) evaluated at /^°^ — > fhcs- 

The time dependence of the reference HCS leads to an explicit time de- 
pendence of L, (, and fhcs- This can be removed by introducing the following 
dimensionlcss variables 



I 



'dt'^, v* = ^-, k*=k£, (30) 
I v (t) 



ff *(k*,v*,*)= fl (k,v,t), C* = -7^, (31) 

where 1=1/ na 3 is proportional to the mean free path and s is proportional 
to the average number of collisions per particle in the time interval (0, t). The 
linearized Boltzmann equation then becomes 
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{O s +ik*-v* +A*)g* = 0, (32) 
where the linear collision operator A* is defined by 

A*Y = L*Y+ C Y /* c -/|- • (vf* hcs Y) , L* = L, (33) 

It is seen that the operator A* is not simply the linearized Boltzmann colli- 
sion operator in dimensionless form, but it includes as well an operator that 
describes the effects of cooling explicitly. 

Normalization of the distribution function requires only that f^ cs g* be 
intcgrable. However, we make the stronger requirement that g* be an element 
of a Hilbert space with scalar product 

(a, b) ee J dv*t hcs (v*)a+(v*)b(v*), (34) 

where a + denotes the complex conjugated of a. The formal solution to the 
linear Boltzmann equation is 

g*(k*,v*,s) = f ^-.e-* s K(z)g*(\i*,v*,0), (35) 

H(z) = (z-ik*-v* - A*)' 1 , (36) 

where the contour encloses the entire spectrum of ik*-v* + A* , both point 
and residual, counterclockwise. All linear excitations are determined from this 
spectrum. The hydrodynamic modes, whenever they exist, are part of the 
point spectrum so it is sufficient to study the eigenvalue problem 

(zk*-v* + A*) ^„(v*) - A n (fc*)Vn(v*). (37) 

The hydrodynamic modes are defined from the macroscopic balance equa- 
tions in the following way. First, a normal form for the cooling rate, pressure 
tensor, and heat flux is assumed to give a closed set of equations (hydrodynam- 
ics). Next, these equations are linearized about the HCS state (n, v = 0, T(t)). 
The deviations of the density, flow velocity, and temperature are then put in 
dimensionless form and given a Fourier representation, 

The resulting linear hydrodynamic equations take the form 

d s Sy a (k, s) + M af 3(k)5y (k, s) = 0. (39) 

The hydrodynamic modes are defined to be the eigenvalues of M a p(k). 
Interestingly, the expression for M a p(k) up through order k does not depend 
on the assumed forms for the pressure tensor and heat flux; only the fact that 
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the cooling rate is normal and dimensional analysis is required. Consequently, 
the corresponding eigenvalues can be considered as known exactly up through 
this order. They are found to be: 

{WW} = jo, £ -£ | + terms of order k\ (40) 



4.1 Existence of Hydrodynamic Modes 

The search for hydrodynamic excitations in the spectrum of the operator 
(ik*-v* + A*) can be performed by assuming they re analytic in k and looking 
first for the k = solutions of (32). The practical issue of constructing these 
modes at finite A: is a separate issue from their existence and is discussed in 
the next subsection. 

Consider again Eq. (21) for the HCS solution in terms of the nonlinear 
Boltzmann collision operator. Although it is an equation for the velocity de- 
pendence of the distribution function, it is parameterized by the hydrody- 
namic fields. Sequential differentiation with respect to these fields introduces 
the linearized Boltzmann collision operator, 

■ [(v u) fhcs] = wJ Mfhcs] = - fhcs{y)L -^r- (41) 

This gives five equations relating the operation of L on derivatives of the 
HCS solution to other derivatives of the same distribution. In detail, these 
terms can be rearranged to construct eigenvalue equations for A* . In this way 
five eigenfunctions and eigenvalues are directly identified [2, 3, 4] 

°.V'-t--T'-t}- < 42) 

«.«»> = {* + V ■ -3 - V ■ 2^., . (43) 

Clearly, these are the hydrodynamic modes in the long wavelength limit de- 
fined by (40). This is therefore a direct confirmation of the existence of hy- 
drodynamic excitations for a granular gas, independent of other derivations 
such as the Chapman-Enskog expansion. 



4.2 Navier-Stokes Order Modes 

Returning to the eigenvalue problem (37) at finite k, the hydrodynamic modes 
can be calculated for small A; by a perturbation expansion 

X n (k)^X^+k\^+k 2 X^ + ---, (44) 

(k) = + + fcV< 2) + ■ ■ ■ , (45) 
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with A„ = A„(0) and ^„ = ip n (0) given by (40) and (43), respectively. 
The analysis is complicated by the three-fold degeneracy associated with the 
unperturbed eigenvalues — (/2. These are the eigenvalues of a vector, which 
can be decomposed into a longitudinal component along k* and two transverse 
components. The operator ik* -v* + A* is invariant under rotations about k* so 
the Hilbcrt space can be decomposed into two invariant subspaces comprised 
of the rotationally invariant functions and their orthogonal complement. The 
transverse and longitudinal eigenfunctions lie in different subspaces. If 
is taken to be the longitudinal component, the eigenvalue problem can be 
separated into two independent eigenvalue problems for n = 1,2,3 and for 
n = 4, 5, respectively. The eigenvalues for the first problem are all distinct, 
and therefore standard perturbation theory applies. The second problem still 
has a two-fold degeneracy that remains even at finite k. 

The fact that ik*-v* + A* is not Hermitian is a second complication. The 

calculations are simplified by introducing a set that is biorthogonal to jvi ^ j, 
W - {l. y + \, k* • v*,e« • v*} , = S nm , (46) 

where jk*,eM;i = 1,2| are three pairwise orthogonal unit vectors. The de- 
termination of the eigenvalues and eigenvectors to order k 2 is then straight- 
forward but lengthy and the details will be given elsewhere [12]. Here, only 
the results for the eigenvalues are quoted. First, all corrections to first order 
in k vanish as expected from the result (37). To order k 2 the coefficients are 

A< 2) = (kv,yl*- 1 k-vV> 1 0) (v)) . 

^=((k.v) a i (A- + ^)" 1 k.^(v)) I 

*S = (V' 2 >-v) (k-v) , (a* + |) _1 k-v^(v^ . (47) 

In the expression of A^ 2 \ Z is a rather involved functional of f^ cs and tp^ 
whose explicit expression will be not given here [3]. It corresponds to the 
second order in the gradients contributions that were neglected upon writ- 
ing expression (17) for the cooling rate in the Navier-Stokes approximation. 
The corresponding results obtained from the linearized Navier-Stokes equa- 
tions are given in terms of the transport coefficients, which are defined there 
from the Chapman-Enskog method. It is possible to show that the results 
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(47) are exactly the same as those from the Chapman-Enskog method, and 
independent confirmation of the latter. 

It has been stressed in section 3 that "hydrodynamics" should be under- 
stood in a more general context than the Navier-Stokes approximation. In the 
present section that means the existence of modes {A„(fc)} beyond quadratic 
order in k. For elastic collisions, it has been shown that the perturbation ex- 
pansion converges [13] so that all higher order terms are meaningful. However, 
the proof does not extend to granular gases. Instead, support for the existence 
of hydrodynamics beyond Navier-Stokes order can be obtained from realistic 
but simpler kinetic model equations. In reference [7] such a kinetic model is 
solved exactly (reduced to quadratures) to demonstrate the hydrodynamic 
modes in the spectrum of the general solution to the linear kinetic equation, 
and to display the equation for the dispersion relations determining {X n (k)}. 
In that case, complete information about the rest of the spectrum is available 
as well, so that the context of a hydrodynamic description can also be studied. 

5 Summary and Outlook 

The origin of hydrodynamics starts with the balance equations. These have 
an exact basis in statistical mechanics and are little more than a statement 
of conservation laws at the macroscopic level. Such equations are correct even 
for very complex fluids, solids, composites, or worse. The next step to obtain 
a closed set of equations for the fields (hydrodynamic, clastic, viscoelastic, 
other) is the point of conceptual and practical difficulty. The choice of con- 
stitutive equations is simple only for systems with well-defined symmetries 
and response characteristics: e.g., simple atomic fluids, elastic solids, liquid 
crystals. In these cases, the symmetries and time scales restrict the form of 
the constitutive equations and the resulting macroscopic equations have little 
uncertainty beyond the numerical values of the constants appearing in them. 

Constitutive equations for complex fluids (e.g., polymers, colloids, foams) 
are more difficult to fix as neither symmetries nor characteristic space and 
time scales are easily determined. Surprisingly, granular gases are more like 
complex fluids than like simple atomic fluids, in spite of their simple structure. 
To a great extent this is due to the single feature of inelastic collisions. Corn 
starch flows like a simple gas under gravity, but responds to large shear as a 
solid. It is probably too ambitious to attempt to find constitutive equations 
encompassing both classes of phenomena. Here attention has been limited to 
conditions of low density gas flow, but even with this restriction it subsumes 
common conditions of both subsonic and supersonic flows for granular gases. 
Thus, it may be expected that simple constitutive relations may have more 
limited applicability than for normal gases. 

In the first part of this presentation, the notion of constitutive equations 
has been tied to the existence of a normal solution to the Boltzmann equa- 
tion. If a normal solution applies, a hydrodynamic description applies (albeit 



14 



James W. Dufty and J. Javier Brey 



with possibly complicated constitutive equations). The idea of a normal so- 
lution is suggested by a two stage process of evolution following preparation 
of the system. The first is a rapid velocity relaxation in each small cell of the 
gas to a distribution characterized by the hydrodynamic fields at each point, 
and also by the gradients of these fields, reflecting the fact that each cell is 
an open system. The second stage is an "equilibration" of the spatial varia- 
tions in the fields from cell to cell, as enforced by the macroscopic balance 
equations. Important open questions at this point include: 1) what is the pre- 
cise definition of a normal solution to the Boltzmann equation, i.e. how is the 
problem properly posed?, 2) is there a "universal" normal solution approached 
for a wide class of preparations?, 3) what is the time scale for approach to 
a normal solution and is it short as compared with the hydrodynamic time 
scales? The Chapman-Enskog perturbation expansion gives partial answers 
to these questions in a very restricted context. Regarding 1), the problem 
is well-posed with solutions assured at each order in the perturbation. The 
Chapman-Enskog normal solution is indeed a universal function of the fields 
and their gradients. Different hydrodynamic states result only from boundary 
and initial conditions determining the values of these fields. Regarding 2) and 
3), the Chapman-Enskog method gives no information. 

The second part of this presentation considers the linear response to small 
spatial perturbations of a homogeneous reference state. The question posed is 
whether the hydrodynamic excitations occur in the spectrum of the generator 
for this linear dynamics. The primary result reported here is that the hydro- 
dynamic excitations do exist in the long wavelength limit (small k). Under 
the assumptions of analyticity in k and convergence of a formal perturba- 
tion expansion, the hydrodynamic spectrum is established. This approach to 
hydrodynamics for a granular gas provides a somewhat sharper mathemat- 
ical formulation of the questions: 1) what is the appropriate function space 
to study the spectrum of (ik*-v* + A*)? and, 2) is the hydrodynamic point 
spectrum found here isolated and smaller than all other parts of the spectrum? 

Simple experiments on driven systems exhibit complex phenomena (e.g., 
symmetry breaking) that nevertheless can be explained by the correct Navier- 
Stokcs hydrodynamics. Then, complex phenomena alone should not preclude 
the simplest Navier-Stokes hydrodynamic description. On the other hand, 
some steady state conditions are necessarily non-Newtonian (beyond Navier- 
Stokes). This should not be interpreted as a failure of hydrodynamics, but 
only as a challenge to explore the correct constitutive equations for these 
conditions. 
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